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Abstract 

The  ^reading  angle  of  an  electrospray  of  Formamide 
with  5-10%  LiCl  by  mass  has  been  measured  to  be  18" 
in  vacuum,  for  conditions  close  to  minimum  stable  flow. 
A  single-stage  einzel  lens  was  able  to  le-focus  the  spray, 
as  predicted  by  paraxial  ray  theory.  In  addition,  a  time- 
dependent  1-D  numerical  model  of  the  cone-jet-spray 
structure  is  described,  and  mitial  results  presented. 

1.  Introduction 

Colloid  thrusters,  working  on  the  principle  of 
electrostatic  extraction  and  acceleration  of  highly 
charged  liquid  droplets,  are  one  of  several  microthruster 
technologies  currently  receiving  attention  for 
applications  ranging  from  main  propulsion  of 
microsatellites  to  high-precision  altitude  control  and 
station  keeping  of  constellation  members.  An 
introduction  to  the  technology  and  a  review  of  its  earlier 
inq)lementation  as  a  mainline  propulsion  concept  (in  the 
1960’s  andl970’s)  can  be  found  in  Ref.  [1]. 

At  MIT,  we  are  currently  pursuing  several  lines  of 
research  on  these  thrusters.  In  this  paper  we  report  on 
our  progress  in  the  direction  of  a  better  understanding 
and  control  of  the  spray  generated  by  a  single  un¬ 
neutralized  colloid  emitter.  One  eventual  goal  is  the 
construction  and  verification  of  models  for  evaluating 
the  possible  interactions  of  these  propulsive  beams  with 
the  spacecraft.  In  addition,  we  aim  at  the  construction  of 
an  improved  Time-of-Flight  mass  spectrometer  for  the 
detailed  study  of  mixed  ion-droplet  sprays,  which  may 
play  a  prominent  role  in  actual  Wgh-performance 
colloid  thrusters  In  order  to  increase  the  resolution  of 
this  device,  ion  optical  concepts  are  needed  for  re¬ 
focusing  of  the  spray  after  its  space-charge  induced 
initial  divergence. 

Section  2  of  this  paper  describes  our  experimental  work 
and  presents  some  data  on  beam  spreading  and 
refocusing  using  a  custom-designed  einzel  lens.  Section 
3  describes  our  on-going  work  on  numerical  modeling 
of  the  colloid  emitter,  starting  with  the  Taylor  cone  and 
following  through  to  the  spreading  electrospray.  Only 
preliminary  results  are  at  present  available,  and  some  of 
these  wiU  be  discussed. 


2.1  Experimental  Apparatus 

The  main  experimental  goal  was  to  measure  the  amount 
of  spreading  m  a  typical  colloid  jet  and  to  determine  if 
the  charged  droplets  can  be  focused  by  means  of 
external  electrostatic  fields. 

The  colloid  emitter  consisted  of  a  single  quartz  needle 
with  ID  of  75fun  and  final  tip  ID  of  30pm.  Fig  1  shows 
a  diagram  of  the  system.  Not  shown  is  the  syringe  pump 
which  was  used  to  inject  and  monitor  the  liquid  flow 
rate  (Formamide  with  5-10%wt.  concentration  of  LiCI) 
into  the  emitter.  Capacitive  effects  in  the  transfer  line, 
however,  did  not  allow  a  precise  determination  of  the 
flow  rate.  Nevertheless  we  know  that  its  value  was 
centered,  within  a  certain  margin,  at  around  2x10'*^ 
mVs. 

The  emitter  is  coated  with  a  metallic  film  that  allows 
electric  contact  between  the  liquid  and  a  high  voltage 
power  supply  (Ve).  The  emitter  tip  was  positioned  at  the 
center  of  a  circular  aperture  of  6  mm  in  diameter  that 
served  as  the  grounded  extractor  electrode.  The 
electrostatic  lens  was  positioned  next  to  the  extractor 
and  was  biased  negatively  with  respect  to  it  (Vp)  by 
using  a  second  high  voltage  power  siqjply. 

The  current  collector  (Fig.  2)  is  comprised  of  a  series  of 
11  concentric  metallic  rings,  each  with  a  width  of  3mm 
(except  the  first  one  which  is  a  circle  with  diameter  of 
1cm)  and  a  maximum  radius  of  3.5  cm.  Each  ring  is 
electrically  isolated  from  the  rest  and  allows  one  to 
measure  current  at  different  radii  with  a  resolution  given 
by  the  number  of  rings.  A  sensitive  amperimeter  was 
used  to  pick-up  the  signal  from  each  ring. 

Experiments  were  performed  inside  a  large  vacuum 
chamber  equipped  with  a  set  of  two  cryogenic  puiq>s  at 
pressures  below  IxlO"^  torr. 

2.2  Beam  Spreading  and  Lens  Design 

One  of  the  most  useful  techniques  to  characterize  the 
beam  of  an  electrospray  is  by  the  use  of  Time-of-Flight 
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Fig.  1 :  Experimental  setiq?  diagram 


(TOP)  spectroscopy  in  which  a  beam  of  emitted 
particles  is  forced  to  drift  towards  an  electrometer  as  the 
extractor  voltage  is  suddenly  cut.  This  allows  one  to 
measure  the  time  it  takes  the  jet  to  reach  the  target  and 
gives  therefore  enough  information  to  estimate  the 
specific  charge  of  the  different  particles  that  compose 
the  beam.  Specific  charge  and  beam  composition  are 
very  important  to  determine  the  performance  of  a 
thruster  as  the  propulsive  efficiency  and  the  specific 
impulse  depend  directly  on  those  parameters. 

Since  the  beam  is  formed  by  charged  particles,  one  can 
expect  it  to  ^read  under  electrostatic  repulsion  forces. 
The  amount  of  this  spreading  will  depend  on  the  beam 
composition  and  its  droamic  characteristics.  From 
previous  experiences  it  is  known  that  spreading 
limits  the  working  distance  for  TOP  to  20-30  cm.  To 
solve  this  problem  one  could  make  use  of  the  well 
developed  theory  of  ion  optics  which  has  been  applied 
for  many  years  to  microfabrication  with  liquid  metal  ion 
sources  ^  \  In  particular,  one  could  make  use  of  a 
focusing  device  called  "einzel”  or  equipotential  lens.  In 
a  lens  of  this  type,  an  accelerating  potential  is  generated 
in  the  path  of  the  charged  beam,  followed  by  a 
decelerating  potential  which  leaves  the  particles  with 
die  same  axial  kinetic  energy  they  had  before  entering 
the  device,  but  provides  a  net  focusing  effect  in  the 
radial  direction. 


Fig.  2:  Current  Collector 


To  calculate  the  beam  trajectory  the  paraxial  ray 
equation  is  used.  In  the  paraxial  approximation  it  is 
assumed  that  the  beam  particles  stay  relatively  close  to 
the  base  trajectory  (beam  axis)  and  ^t  the  beam  profile 


has  a  uniform  charge  distribution.  The  force  equations 
can  be  solved  considering  axially  symmetric 
electrostatic  fields  with  cylindrical  geometry.  These 
fields  are  expressed  as  the  first  order  terms  of  the 
solution  in  series  of  Poisson's  equation.  What  is  most 
interesting  about  the  paraxial  approximation  is  that  we 
only  require  knowledge  of  the  electrostatic  potential 
{Vo)  at  the  axis  of  symmetry.  The  paraxial  ray  equation 
can  then  be  written  as 


^  Vy  Por 

^  o  I  o  —  ___  *  o 


2V  4V 


(1) 


where  the  space  charge  density  po  is  given  for  'N  species 
in  the  beam  as 

/i  and  {q/rn)i  are  the  current  and  the  charge  to  mass  ratio 
of  the  rth  species  in  the  beam.  (A  single  species  was 
assumed  in  this  paper).  The  on-axis  potential  can  be 
extracted  after  solving  L^lace's  equation  for  a  given 
geometry  of  electrodes  with  the  help  of  an  elliptic 
solver. 


Fig.  3  shows  the  geometry  of  the  emitter,  extractor  and 
the  lens,  superimposed  to  the  mesh  used  to  solve 
Laplace's  equation,  and  the  solution  itself.  It  is  clear 
bow  the  strong  negative  potential  m  the  middle  region 
of  the  lens  is  effectively  shielded  out  by  the  grounded 
electrodes  positioned  in  the  inlet  and  outlet  planes.  The 
next  step  is  to  extract  the  axial  potential  from  this 
configuration  and  use  it  to  solve  (1)  with  the  conditions 
expected  during  the  focusing  experiment.  Fig.  4  shows 
the  result  of  the  calculation.  This  particular  beam  shq>e 
was  obtained  with  a  negative  potential  Vp  =  5.5  kV.  The 
rest  of  the  required  parameters  were  extracted  from  the 
experimental  conditions  to  be  presented  in  the  next 
subsection. 


23  Experimental  Data 

The  first  series  of  experiments  were  aimed  at  estimating 
the  amount  of  spreading  of  the  charged  beam.  The 
results  are  summarized  in  Fig.  5,  where  current  density 
measurements  were  performed  at  three  different  axial 
positions  from  the  extractor/emitter  plane.  It  is 


Fig.  3:  Computational  domain  -  emitter,  extractcM*  and  einzel  lens 


Fig.  4:  Trajectory  calculation  from  the  paraxial  ray  equation 


interestmg  to  note  fliat  for  the  first  two  measurements  (L 
=  3.5  and  7  cm)  a  more  or  less  xmiform  charge  density 
profile  can  be  observed.  For  the  third  measurement  (L  = 
13  cm)  the  profile  is  also  almost  uniform  but  occiqjies 
the  whole  area  of  the  detector.  From  the  cutoff  values  of 
the  current  and  the  distances  involved,  a  semi-angle  of 
about  18  *  can  be  calculated. 

In  the  second  experiment,  the  einzel  lens  was  installed 
in  the  system  and  the  current  collector  was  positioned  at 
a  distance  L=13  cm.  Fig  6  shows  the  current  collected 
by  the  innermost  circle  of  the  detector  (r  <  0.5  cm)  as  a 
function  of  various  negative  lens  potentials.  We  observe 
that  a  maximum  in  current  (72  nA)  is  obtained  when  Vp 
is  around  5.5  kV.  This  value  of  current  was  introduced 
in  the  solution  of  die  paraxial  ray  equation  in  subsection 
2.2.  The  result  of  the  calculation  indicates  that  a  beam 
"spot"  of  less  than  0.5  cm  in  radius  is  generated  at  the 
distance  where  the  detector  is  (Fig.  4.) 

In  all  cases  the  emitter  to  extractor  potential  was  held  at 
1.9  kV.  This  value  was  determined  at  atmospheric 
pressure  by  visual  inspection  of  a  stable  cone  at  the 
needle  tip  and  measurement  of  a  stable  current. 

2.4  Discussion 

Experimental  results  have  indicated  that  focusing  is 
possible  for  beams  of  charged  particles  emitted  firom 
colloid  thrusters.  Nevertheless,  more  work  is  required  to 


refine  the  results  presented  here.  In  particular, 
controllability  of  the  flow  rate  is  important  for  both  the 
determination  of  some  spray  characteristics  and  for 
ensuring  good  experimental  repeatability. 

Even  Avith  such  uncertainties,  the  measurements  are  well 
correlated  to  mathematical  modeling.  Future  work  will 
include  a  refinement  in  the  measurements  presented 
here,  including  TOF  spectroscopy  with  the  help  of 
electrostatic  focusing. 


Fig.  5:  Current  density  distributions  (experimental) 


Enz^-tens  charactenalics 


Fig.  6:  Electrostatic  focusing  (e;q>enniental) 


3.  A  Dynamic  Model  cf  an  Emitting  C6ne-Jet 

The  structure  and  performance  of  the  individual 
emitting  cone-jets  which  underlay  colloid  thrusters  (and, 
in  sli^tly  modified  form,  field-emission  thrusters)  has 
been  the  object  of  many  studies  by  electro-spray  and  ion 
enussion  spectroscopy  scientists.  The  seminal  work  was 
done  by  Taylor  ^o  first  identified  the  49.3°  cone  as 
die  electro-capillary  equilibrium  shape  under  idealized 
conditions.  The  equilibrium  shapes  and  internal 
pressures  for  more  general  conditions  were  conqiuted 
by  Pantano  and  Gahan  Calvo  who  used  a  boundary 
element  method  for  the  (axisymmetric)  electrostatic 
calculations.  They  imposed  Taylor’s  angle  as  a  tip 
boundary  condition  (in  place  of  the  thin  jet  which 
actually  develops),  and  were  able  to  determine  the 
cone’s  volume  and  existence  voltage  as  a  function  of 
size  and  internal  pressure.  Gahan  Calvo  ^  developed  an 
analytical  model  for  the  cone-jet  combination  with  a 
few  reasonable  simplifications,  and  obtained  scaling 
laws  and  definite  predictions  for  jet  ^ape  and  diameter, 
surface  field  distributions  and  emitted  current.  Similar, 
more  detailed  analytical  results  were  presented  by 
Chemey  ^  recently.  Both  Gahan  Calvo ’s  and  Chemey’s 
results  confirm  and  explain  the  semi-empirical  laws 
obtained  by  Fernandez  de  la  Mora  based  on  extensive 
experimentation  with  many  conducting  fluids. 

Our  current  work  is  aim^  at  extending  these  previous 
results  by  constructing  a  flexible  numerical  model  of  the 
cone-jet  combination,  which  can  serve  as  the  basis  for 
assessing  a  large  number  of  parameters  and  physical 
effects.  The  formulation  is  axisymmetric,  and  quasi-one 
dimensional,  but  not  necessarily  steacfy;  solutions  are 
sought  by  time-marching,  which  is  helpful  for 
understanding  the  numerical  issues,  and  also  yields 
information  on  the  dynamics  of  the  cone-jet  structure. 
Much  of  the  formulation  was  initiated  by  V.  Khayms 
who  pursued  a  steady  state  solution  through  direct 
iteration. 


3.1  Formnlation 

Fig.  7  ^ows  a  schematic  of  a  section  of  the  liquid.  The 
volume  in  an  element  such  as  that  shaded  in  the  figure, 
divided  by  dx,  is 


Vol- 


cosa(l  +  cosa) 
and  the  volume  flux  crossing  it  is 


Q.= 


ItcR^u 


1  +  cosa 

The  mass  conservation  law  is  then 


dm  ,  dQ, 

dt  dx 


(3) 


(4) 


(5) 


Fig.  7;  Modeled  Geometry.  The  shaded  volume  element  is  a 
section  cf  a  q^iencal  shell  perpendicular  to  the  liquid  sudTace.  AH 
^JJanthies  are  assumed  constant  on  one  such  shell,  and  u,  £t  are 
assumed  perpendicular  to  its  ^herical  boundaries. 

Similarly,  the  net  charge  in  the  element  (per  unit  axial 
length),  is 


Ch  = 


cosa 


(6) 


wliere,  if  the  fluid  is  conductive,  a,  =  eJSny  E„  being  the 
outside  normal  electric  field.  The  bulk  carries  no  net 
charge.  The  net  charge  flux,  or  current,  contains  two 
contributions,  one  from  surface  charge  convection,  the 
other  from  bulk  conduction.  The  fluid  conductivity,  Ky 
is  treated  as  a  constant,  and  one  obtains 


V  1  +  cosa  J 


(7) 


Hraice,  charge  conservaticm  requires 


d{Ch)  dl, 

dtdx 


\^ere  possible  ion  emission  or  electron  capture  at  the 
surface  have  been  ignored.  The  x-momentum  per  unit 
axial  distance  is 


{Mom)  = 


puTtR^ 
cos  a 


(9) 


and,  including  the  electrical  surface  traction  Xe  - 
die  axial  momentum  balance  is  then 

d{Mom)  ^  ^  ^  iTtRBQE^E^  (10) 

at  dx  dx 


The  final  mechanical  equation  required  is  a  balance  of 
normal  forces  on  a  surface  element: 


^0  r>2  cosa 

^  2  ”  R 

/*  dR^ 

(tana  =  -— ) 
dz 


3  d^R 

-cos  ce- 


dj^ 


(11) 


The  full  matrix  is  then  partitioned  into  blocks 

"  representing  potential  on  liquid  points  due  to  charges  on 
liquid  elements  {N^^X  potential  at  liquid  points  due  to 

charges  on  the  metal  supply  tube  {N^^),  etc.  We  dien 
have  the  matrix  equations 

(P'=N,,E/  +  N„E:  (16) 

(j>‘=N„E,^  +  N„Ej  (17) 

^ere  are  numerical  vectors  collecting  all 

normal  fields  on  the  licpiid  and  on  die  tube,  respectively. 
Inqiosing  a  known  tube  potential  K,  and  eliminating 
be^een  (16)  and  (17)  leads  to 

<l>‘  =  N,,N-^VHNu-N„N„-%,)EJ  (18) 

'vshich  is  the  desired  (global)  relationship  between 
normal  fields  and  potentials  on  the  liquid  surface.  The 
numerically  intensive  matrix  inversion  needs  to  be 

done  only  once  for  a  given  geometiy,  but  the  matrices 
have  to  be  iqidated  as  the  liquid  shape 

evolves. 


where  die  two  terms  on  the  right  are  the  two  principal 
curvatures  of  die  liquid  surface. 

The  tangential  field  is  related  to  the  surface  potential 
throu^ 

E=-cosa^  (12) 

dx 

The  relationship  between  0  and  the  normal  field 
E^aje^  is  more  complex,  since  every  charge  in  the 
^stem  (on  the  liquid  as  well  as  on  the  electrodes) 
influences  each  local  potential.  We  use  for  this  a 
boundary-element  method  similar  to  that  in  Ref.[6].  The 
contribution  S^t  to  the  potential  at  point  f,  due  to  the 


diarge  on  a  surface  ring  at  J,  with  axial  extoit  Sx,  is 

written  as 

II 

(13) 

where 

”  2fr  '  cosa,. 

(14) 

and 

n, 

(15) 

and  K(m)  is  the  complete  elliptic  integral  of  the  first 
kind  hi  order  to  simulate  a  grounded  electrode  at  x  =  Z, 

an  “image  ring”  with  an  equal  and  opposite  charge  is 
also  included.  This  amounts  to  replacing  '\fmK{m)  by 

.^K{mi)  in  (14),  ■wdiere 
but  with  Xj  replaced  ly  2L—Xj  . 

is  as  in  (15), 

ZJ.  Numerical  Implementation 

A  fixed  1-D  grid  is  generated,  with  highest  resolution 
near  the  expected  cone  tip  and  near  the  tube  mouth, 
where  a  field  singularity  is  expected  due  to  the  sharp 
change  of  surface  slope.  For  initialization  of  variables, 
as  well  as  for  selection  of  appropriate  grid  spacing,  the 
dimensionless  results  of  Ref  [7]  have  been  used.  Spatial 
derivatives  are  approximated  to  order  on  the  uneven 
grid,  and  a  first-order  Euler  st^  is  used  for  tenqmral 
updating  of  (Fb/),  {Ch)  and  {Mom\  using  Eqs.  (5),  (8) 
and  (10),  respectively.  After  the  updates,  the  primary 
variables  (/?,a,ii,£^„)  are  extracted  from  the  iq>dated 
quantities,  p  is  calculated  from  (11),  and  £q.  (18)  is 
used  to  calculate  the  new  potential  distribution.  E^  then 

follows  from  (12),  and  a  new  updating  cycle  can  begin. 

The  time  step  Af  used  for  iqxlating  needs  to  be  smaller 
than  both,  the  charge  relaxation  time  =  e^lK  and 

the  smallest  passage  time  through  one  grid  cell, 
~  (^1  liquids  of  high  conductivity,  {K  > 

1  Si/m),  is  less  than  about  10”^^  s,  well  below 
and  short  enough  that  complete  simulation  through  a 
few  microseconds  may  be  impractical.  For  these 
situations,  a  modified  procedure  is  necessary,  in  which 
the  cone  region  is  assumed  electrically  (not 
mechanically)  relaxed  all  the  time  (i.e,  /^(jc,0  =  /(0)> 

but  the  jet  region  current  is  allowed  to  evolve 
(fynamically.  The  level  of  current  I(t)  is  then  determined 
hy  the  value  of  at  the  cone-jet  transition.  This  has  not 

yet  been  implemented  Instead,  work  has  concentrated 


on  cases  with  lower  conductivity  (ir<0.01  Si/m),  for 
which  an  extensive  data  base  exists  Refs.  [2^]. 


3.4  Jet  Instability  and  Spray  Dynamics 


Even  in  the  absence  of  charge,  the  jet  would  pinch  itself 
into  droplets  at  some  downstream  distance,  following 
growth  of  the  Rayleigh-Taylor  instability.  The  presence 
of  charge  can  only  accentuate  this  tendency.  For  purely 
convected  surface  charges,  a  linearized  stability  analysis 
presented  in  Ref.[10]  obtained  die  dispersirm  relation 


(O' 


^ _ r  x/,(x)r  1 

1,{X)  \  1  /,(jt) 


(19) 


^ere  /is  the  surface  tension,  is  the  local  mean 
radius,  x  —  kR^ ,  k  being  the  wavenumber,  and  /j  are 
modified  Bessel  functions,  and 


is  a  measure  of  the  state  of  charge,  \^ich  is  of  order 

unity  in  electrospray  jets.  Negative  values  of  (O  imply 
die  existence  of  an  exponential  growth  mode.  Eq.  (19)  is 
plotted  in  Fig.  8.  The  case  =  0  is  the  classical  result 

of  Raylei^  and  Taylor,  with  peak  growth  at  =  0.7 , 
leading  to  an  eventual  radius  of  the  pinched  droplets 
Rjy  =  For  FV  ~  1,  somei^at  smaller  droplets 

are  indicated  (i?^  ^1,6 -1.7/?;^),  but  also  faster  growth 
of  the  instability.  Also  of  interest  for  our  simulation,  all 
wavelengths  are  seen  to  be  unstable,  although  only 
weakfy  so  for  the  longer  ones  (small  x). 


In  the  jet  simulation,  the  axial  grid  spacing  Ax  is  of  the 
order  of  the  jet  radius  at  the  neck,  ahhou^  it  is  later 
allowed  to  grow  much  larger.  It  is  dierefore  possible  to 
capture  the  growth  of  the  jet  instability  and  the  jet 
pinching  near  the  neck.  Fig.  9  shows  an  example, 
corresponding  to  a  Formamide  cone-jet  with  AH).  01 


Si/m,  2  =  7.84x10"*^  mVs,  F=2000V.  The  time  after 
simulation  start  is  t  =  5AxlO"’s,  and  x  =  1.2xl0”^m 
corresponds  to  about  3-4  jet  diameters  downstream  of 
the  cone’s  neck. 


Using  y=  0.05N/m,  p  =  113(Scg/m^  and  (from  Fig.  8)  a 
maximum  normalized  (-m^)  of  0.5,  the  exponential 
time  constant  for  growth  is  predicted  to  be  7  x  10“*^, 
and  the  most  unstable  wavelength  to  be  3pm.  Fig.  9 
shows  pinching  in  about  8  time  constants  (a  reasonable 
result),  and  rou^ly  conftrms  the  3pm  wavelength. 

A  detailed  computation  of  the  jet  breakip  process  is 
however,  beyond  the  scope  of  this  quasi-one-D  model. 
In  order  to  continue  capturing  the  electrostatic  effect  of 
die  jet  on  the  cone,  we  need  to  replace  the  jet  equations 
by  a  set  of  “spray”  equations  beyond  the  breakup  point 
For  this,  a  time-dependent  variation  of  the  Paraxial  Ray 
method  will  be  used,  although  this  is  yet  to  be 
inplemented.  The  spray  flux  is  assumed  uniform  xsp  its 
edge  (consistent  with  our  eiqperimental  results  else^ere 
in  this  paper),  and  the  expansion  yields  a  radial  field  at 
r^R  of 


E  M 

”  2  dx^  (qlm)2jcRe, 


where  M  is  the  mass  contained  in  unit  length  of  the 
spray  (and  it  obeys  a  conservation  equation 

=  0),  and  q/m  is  the  locked-in  charge-to- 

dt  dx 

mass  ratio,  carried  over  from  the  jet.  The  velocity 
components  n(axial)  and  v(radial)  obey 


du  du 
dx: 


—E, 


4  +  = 

dt  dx  m 


(22) 

(23) 


and  the  spray  radius  then  evolves  according  to 


dR 

dt 
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(24) 
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Rg.  10:  Preliminaiy  numerical  results.  Dotted  line,  r  =  0. 
Continuous  line,  r  =  1.3  jjis 


With  some  uncertainty,  arising  from  the  details  of  the 
breakiq)  process,  tiiese  equations  will  allow  conq>utati(m 
of  the  spread  of  non-neutralized  electrospray  plumes, 
and  will  find  use  in  the  design  and  analysis  of  TOF 
e7q>eriments,  as  planned  in  our  laboratory. 

3^  Some  Preliminary  Results 

Fig.  10  shows  a  sample  of  results  pertaining  to  a  cone- 
jet  of  Formamide,  with  a  conductivity  XK).  001  Si/m,  at 
a  flow  rate  (^7.84x10*^*  mVs,  which  is  a  few  times  the 
minimum  stable  flow  for  that  conductivity  ^"*1  This 
conductivity  level  is  also  what  can  be  obtained  in 
concentrated  Glycerol  solutions,  as  in  the  first 
generation  of  colloid  thrusters  The  delivery  tube,  at 
F=2000Volt,  has  a  radius  of  10  pm  with  a  sharp, 
chamfered  1^,  and  is  assumed  placed  Imm  away  from  a 
grounded  collector  plane.  Shown  are  profiles  of  liquid 


radius  /?,  velocity  w,  local  current  /„  pressure  p,  normal 
electric  field  and  potential  For  each  of  these,  we 
show  the  initial  profiles  assumed,  and  those  conqiuted 
after  about  1.3xl()“®5.  This  time  is  about  1/8  what  it 
would  take  to  "fill"  a  Taylor  cone  of  die  assumed  size, 
at  the  assumed  flow  rate,  so  that  convergence  is  still 
inconqilete.  The  time  step  used  was  2xl0'*‘5,  about 
1/500  of  the  fluid’s  relaxation  time;  increasing  this  to 
2  x  10"*®5  resulted  in  strong  numerical  oscillations  near 
the  tube  exit 

Several  observation  can  be  made:  Firstly,  at  1.3  ps,  the 
local  current  /,  has  relaxed  to  a  near-constant  value  of 
8x10“^  A,  several  times  higher  than  the  expected 
steady-state  value,  based  on  the  empirical  data  of  [4], 
which  was  used  in  the  initialization  (and  is  visible  in  the 
parts  of  the  not  yet  affected).  This  level  exhibits  a 
slow  decay  rate  over  time,  as  the  extent  of  the  relaxed 
zone  spreads  forward.  It  should  be  noted  that  the 


flattening  of  the  Ij  profile  was  already  evident  at  times 
of  the  order  of  1x10"*^,  i.e.,  the  relaxation  time  at  the 
diosen  conductivity. 

The  peak  normal  field  is  about  0. 2-0.3  V/nm^  and  occurs 
near  the  front  of  the  spreading  relaxation  zone.  Because 
of  its  hi^  value,  the  liquid  pressure  is  seen  to  change 
from  its  initial  distribution  (near  zero  in  the  cone, 
positive  in  the  jet)  to  between  1  and  2  atm.  negative. 
This  might  inq)ly  liquid  breakiq),  althou^  metastability 
may  hold  it  togedier  during  this  transient.  The  potential 
is  nearly  constant  in  the  cone,  falling  by  some  200  Volt 
near  the  tip,  and  more  rapidly  beyond  it  The  cone  itself 
has  been  elongating,  from  some  9  pm  initially 
(following  the  assumed  Taylor  angle),  to  some  20  pm  at 
die  1.3  ps  time  point 

We  note  finally  that  the  jet  has  not  yet  ^own  signs  of 
pindiing  in  this  simulation,  whereas  for  ^=0.01  Si/m, 
pinching  was  ^own  (Fig.  9)  at  about  0.5  ps.  This  is  to 
be  expected  from  the  dependence  of  the  instability 
growdi  time  (Eq.  19),  and  the  variation  of  at  minimum 
flow  as  K"^  (Rrf.  4). 

4.  Conctnskiiis 

The  ability  to  reverse  the  electrostatic  spreading  of  a 
non-neutralized  colloidal  jet  has  been  experimentally 
demonstrated,  which  opens  the  way  for  improved 
diagnostic  instrumentation.  A  new  method  for 
numerical  simulation  of  the  cone-jet  dynamics  has  been 
developed,  and  preliminary  results  are  encouraging. 


universal  scaling  laws  in  electrospraying^  Phys. 
Rev.  Let.,  V.  79,  N.  2,  july  1997. 

8.  L.  Qhemey,  “Structure  of  Taylor  cone-jets:  limit  of 
low  flow  rates’’.  J,  of  Fluid  Mech.,  (1999).  V.  378, 
pp  167-196. 

9.  V.  Khayms,  "Advanced  Propulsion  for  Micro¬ 
satellites",  Doctoral  Thesis,  MIT,  Sept  2000. 

10.  M.  Martinez-Sanchez,  et.al.,  in  the  Mterim  Progress 
Report  on  NASA  SBIR  16.01,  May  2000. 


1.  M.  Martinez-Sanchez,  J.  Fernandez  de  la  Mora,  V. 

Ifruby,  M.  Gamero-Castano  and  V.  Khayms, 
'Research  on  Colloid  Thrusters”.  Paper  lEPC  99- 
014,  26^  International  Electric  Propulsion 

Ccmference,  Kitakyu^u,  Japan,  Oct.  1999. 

2.  M.  Gamero  and  V.  Hruby,  "Electrospray  as  a  source 
of  nanoparticles  for  efficient  colloid  thrusters", 
paper  AIAA-2000-3265,  36*  Joint  Propulsion 
Conference,  HuntsviUe,  AL,  July  2000. 

3.  J.W.  Ward,  R.L.  Seigler,  "Trajectory  calculations  of 
liquid  metal  ion  sources",  J.  Vac.  Sci.  Technol., 
19(4),  Nov/Dec  1981,  p.  1082 

4.  J.  Fernandez  de  la  Mora  and  J.G.  Loscertales,  “The 
current  emitted  by  highly  conducting  Taylor  cones”. 
J.  Fluid  Mech.,  vol  260,  pp  155-184  (1994) 

5.  G.I  Taylor,  'TMsintegration  of  water  drops  in  an 
electric  field”,  Proc.  R.  Soc.  London,  A  280.  383- 
397. 

6.  C.  Pantano,  AM.  Gahan-Calvo  and  A.  Barrero, 
“Zeroth-order  electrostatic  solution  for 
electrospraying  in  Cone- Jet  mode”  J.  Aerosol  Sc.,  V. 
25,N.  6, 1065-1077  (1994). 

7.  A  Ganan-Calvo,  “Cone- Jet  analytical  extension  of 
Taylor’s  electrostatic  solution  and  the  asymptotic 


